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Starting from a Hamiltonian description of the photon within the set of Bargmann-Wigner equa- 
tions we derive new semiclassical equations of motion for the photon propagating in static gravita- 
tional field. These equations which are obtained in the representation diagonalizing the Hamiltonian 
at the order h, present the first order corrections to the geometrical optics. The photon Hamiltonian 
shows a new kind of helicity-torsion coupling. However, even for a torsionless space-time, photons 
do not follow the usual null geodesic as a consequence of an anomalous velocity term. This term 
is responsible for the gravitational birefringence phenomenon: photons with distinct helicity follow 
different geodesies in a static gravitational field. 



PACS numbers: 



I. INTRODUCTION 



In the last few years many studies focused on the trans- 
port of quantum particles with spin. Indeed, manipulat- 
ing spin polarization of electrons is a challenging goal in 
semiconductor spintronics. Achieving this goal requires 
the understanding of the spin transport mechanism in 
systems with spin orbit (SO) interaction. It was found 
that in such a system a Berry phase in momentum space 
plays an important role by affecting both particle phase 
and its transport properties It is well known since 
the seminal work of Berry [2j , that when a quantum me- 
chanical system has an adiabatic evolution, a wave func- 
tion acquires a geometric phase. It is only recently that 
the possible influence of the Berry phase on transport 
properties (in particular on the semiclassical dynamics) 
of several physical systems has been investigated. Semi- 
conductors, having SO couplings greatly enhanced with 
respect to the vacuum case, require a theory of spin trans- 
port. However, even in the vacuum, new fundamental 
results concerning the semiclassical equations of motion 
of electrons were recently derived. For instance in [3| and 
[H, considering the Dirac equation in an external poten- 
tial, it was shown that the position operator acquires a 
spin-orbit contribution which turns out to be a Berry 
connection rendering the coordinate algebraic structure 
non-commutative. This drastically modifies the semiclas- 
sical equations of motion and implies a topological spin 
transport similar to the intrinsic spin Hall effect in semi- 
conductors A similar non-commutative algebra has 
been also found in the context of electrons in magnetic 
Bloch bands [5| , leading to an anomalous velocity term. 

Despite its very different nature, the photon displays 
many similar behaviors with electronic phenomena such 
as energy bands in photonic crystals and localization. 
These similarities stem from the wave-like nature of 
quantum particles. Because a photon is also a spinning 
particle it is important to understand if SO interaction 
may influence the transport of light in a similar way as 
electrons in vacuum or in semiconductors. It has been 



known for long that there is no position operator with 
commuting components and as a consequence photons 
are not localizable. Therefore, one of the main differ- 
ences between photons and electrons also disappears as 
the coordinates of both particles have in fact noncommut- 
ing components. It is precisely this property which is at 
the origin of the SO interaction leaving no doubt about 
its contribution to the propagation of photons. There- 
fore these recent studies of the SO coupling in different 
systems taught us that one can treat both kind of parti- 
cles on an equal footing. It is thus even more legitimate 
to wonder whether electronic phenomena have photonic 
counterparts. 

In addition, the localization of light rays is the essential 
ingredient of the construction of Minkowski space-time. 
Also in the context of general relativity, it has been ar- 
gued that the Riemannian metric is determined by the 
properties of light propagation [6] 7] . A deeper under- 
standing of the properties of light and in particular its 
SO interaction which is at the origin of the non localis- 
ability, is thus necessary to build a quantum version of 
space-time. 

It has already been observed that the SO coupling in- 
duces a rotation of the polarization plane of light prop- 
agating in an optical fiber with torsion This effect 
had been predicted long ago by Rytov and Vladimirskii 
3] and can be interpreted in terms of Berry phase 
2| in momentum space. Recently it was shown that the 
noncommutativity of the coordinates affects the ray of 
light itself in an isotropic inhomogeneous medium [H,[llj]. 
Other approaches to the spin-orbit contributions for the 
propagation of light in isotropic inhomogeneous media 
have been the focus of several other works [12| and have 
led to a generalization of geometric optics called geomet- 
ric spinoptics [HI]. 

In this letter we investigate how the photon propaga- 
tion in a static gravitational field, which can be seen as 
an anisotropic inhomogeneous medium, is affected by its 
spin-orbit interaction. Our photon description is based 
on a Dirac like Hamiltonian in an arbitrary static grav- 
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itational field with a possible (non-specified) torsion of 
space. The Bargmann-Wigner equations allow us to 
build the wave function of a spin one particle and to 
deduce the dynamical operators which satisfy unusual 
commutation relations in the representation where the 
Hamiltonian has been diagonalized at the order Ti. It 
is in the diagonal representation (similarly to the Foldy 
Wouthuysen (FW) representation) that the physical con- 
tent of the theory is best revealed. In particular we find 
that the photon semiclassical Hamiltonian shows a kind 
of helicity-torsion coupling resulting from the presence of 
Berry curvatures. From the semi-classical Hamiltonian 
we can deduce new equations of motion taking Berry- 
curvatures into account. These equations correspond 
to the first corrections to the geometrical optics. It is 
found that the helicity is always conserved and that even 
in a torsionless space-time, photons do not follow the 
null geodesic due to an anomalous velocity term. The 
presence of this term is a general feature of the modern 
point of view of the topological spin transport of quan- 
tum particles [l[ Q [J] [lj| . The anomalous velocity is di- 
rectly responsible for the gravitational birefringence phe- 
nomenon, i.e., photons with distinct helicity follow differ- 
ent geodesies in a static gravitational field. This confirms 
earlier claims suggesting that the absence of birefringence 
in Einstein's gravity might be only a consequence of the 
limitation of the geometrical optics limit [151 ]. Our com- 
putation also shows that the velocity of light is still equal 
to c at this order of approximation (a polarization depen- 
dent velocity would be a T% 2 phenomenon). Finally, it is 
worth noticing that there exists another spinorial repre- 
sentation of the photon Hamiltonian equivalent to the 
Maxwell equation [T(| which contrary to the Bargmann- 
Wigner equations does not suffer inconsistencies due to 
the redundancy of the photon wave functions. However it 
can be seen that both approaches lead to the same semi- 
classical energy and dynamical operators. The advantage 
of the Bargmann-Wigner approach presented here is that 
it is more usual, more straightforward and that it removes 
also the null energy eigenstates. 



irreducible representation of angular momenta 1 deduced 
from the composition of two states with angular momen- 
tum 1/2. It can be proved that Eqs[T] are equivalent to 
the Maxwell equations [Isj . 

Next, we write an Hamiltonian associated to each 
Bargmann-Wigner equations = a^KP which can be 
diagonalized by the product of usual FW unitary trans- 
formations [lg, [/W(P) = (E + c(3a^ r \P)/EV2 such 
that 

^(i)(P)^(2)(p)#M [/ (i)(p)+ [/ (2)( P )+ = E p(r) t ( 2 ) 

with E = Pc the energy of the photon in the vac- 
uum. The wave function is also transformed and be- 
comes * (aioa) (P) = E- 1 (U^(P) + U^(P) + )^ {aia2) (P) 
in the FW representation. Whereas the quasi- 
momentum is invariant through the action of U^ r \ i.e., 
p=[/( 1 )j7( 2 )pZ7( 1 )+£/'( 2 )+ = p ; the position operator be- 
comes r = ihdp + A'-- 1 ' 1 + A^ with a pure gauge po- 
tential A {r) = ihU {r \p)d p U {r \p)+ induced by the FW 
transformation. After projection of r on the positive en- 
ergy subspace we get the non trivial SU(2) connection 
A =(p A S)/p 2 , and a new position operator r for the 
free photon r = R + (p A S)/p 2 where R =ihd p is the 
canonical coordinate operator and £ = %(a^ + a^ )/2 
the spin one matrix. This definition of the position op- 
erator gives rise, through the commutation relations, to 
a monopole in momentum space [a; 1 ,^] = iTi6 lJ (p) = 
—ihe^ k Xpk/p i where A = ±7i is the helicity and (p) = 
V p iAj — VpjAi + [Ai t Aj] the so called Berry curvature 
(the origin of this name is explained in see also [20|). 

We then build the Hamiltonian of the free photon 
of positive energy as H = (J3W + E^)/2 = pc with 
£?( r ) = pc the positive eigenvalue of the operator H^ r \ 
Therefore, for free particles, the dynamical equations of 
motion in the FW representation are trivially given by 
dv/dt = &[t,H] = pc/p and dp/dt = j-[p,H] = which 
in particular implies for the light velocity v — c. Note 
that the same results can be deduced from the Maxwell 
equations [2l| . 



II. FREE PHOTON 

We start with the description of the photon dynamics 
in the vacuum by considering the two Bargmann-Wigner 
equations (r = 1, 2) of a massless spin 1 particle [l7j |: 

(l i o ) d o + 1 ( r) d l )y (aia2) (x)=0 (i= 1,2,3) (1) 

where>I'( aia2 ) = ^f a2 ai) 1S the symmetrized Bargmann- 
Wigner amplitude with at running from 1 to 4 the spino- 
rial indices of the wave function (as usual the indices 
and i refer respectively to the local time and space co- 
ordinates, with the flat metric (+,—,—,—) ). The 7*- r ' 
matrix is a 7 matrix acting on a\. This symmetrized di- 
rect product of two Dirac spinors assures that the positive 
energy subspace forms a 3-d space corresponding to the 



III. PHOTON IN A STATIC GRAVITATIONAL 
FIELD. 

We now extend our previous approach to the case 
of a photon propagating in a arbitrary static gravita- 
tional field, where <?oi = for i = 1,2,3, so that 
ds 2 = goo(dx ) 2 — gijdx l dx : > = 0. Consider again the 
Bargmann-Wigner equations of motion with the follow- 
ing associated Hamiltonian of the Dirac form 

H {r) = V5oW r) -P + ^rgV + *~r°% (3) 

and P given by P ct =/i 4 Q (R)(P J + |e e/37 rf ai) with h\ the 
static orthonormal dreibein (a = 1,2,3), T" 13 the spin 
connection components and e Q /3 7 cr 7 = |(7 Q 7 /3 — 7 /3 7 Q ). 
The coordinate operator is again given by R =ihd p . Note 
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that here we consider the general case where an arbitrary 
static torsion of space is allowed. It is known [12] that 
for a static gravitational field (which is the case consid- 
ered here), the Hamiltonian is Hermitian. We now 
want to diagonalize H ^ through a unitary transforma- 
tion U^(P). Because the components of P depend both 
on operators P and R the diagonalization at order h is 
performed by adapting the method detailed in [23[ to 
block-diagonal Hamiltonians. To do so, we first write 
in a symmetrical way in P and R at first order in 
h. This is easily achieved using the Hermiticity of the 
Hamiltonian which yields 
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Now for a ease of exposition we temporarily include the 
goo is the definition of the vierbein h % a — > gooh % a . The 
contributions of the goo will be again explicitly written 
later on in the final result for the energy. 



A. Semiclassical Hamiltonian Diagonalization 

The semiclassical diagonalization is achieved in two 
steps. In the first one we diagonalize at first order in 
Ti by considering the formal situation where R is consid- 
ered as a parameter commuting with P. In the second 
one we add the contributions due to the non commuting 
character of the dynamical variables. 



1. Diagonalization when P and R commute. 

The Hamiltonian Hq (we add the index when R is 
a parameter) can then be diagonalized at first order in h 
by the following unitary FW matrix 



U^(P) = 



>/2«) 2 V 2 V 

+P + .a (r) l +N) 



(4) 



with E { r) = 



P+P .qMg 



We introduced 



also the notations N = f - -^ X1 ^ and D = 1 + 

f ^ ^ Px 2 ^?3 xP with r 07 = e e p^V^. (here, and in the 
i 

sequel, P means (P 2 ) 2 , with the metric and P 3 is 

3 

obviously defined as (P 2 ) 2 )■ 

The proof of this diagonalization relies on the following 
properties: For each parameter R the matrices h % a and 
r"' 3 are independent of both the momentum and position 
operators. The matrices and a.P anticommute and in 

(r) 

the Taylor expansion of Eg all terms commute with p 
and a^.P + P + .a^. In this context the diagonalized 



Hamiltonian is equal to H^U^ + which reads 



U {r) H {r) U [r)+ 
u n n o u o 



cB 



&(PX„)P.<7< r > 



4 p2 



(•5) 



which can also be written 



U^H^U^ = c^\jP^+'±P.a^e e01 Tf^ 



+ 



s(P.r )P.aW 



4 p2 



(6) 



In this last expression we have neglected contributions of 
curvature type of order Ti 2 . 



2. Corrections when P and R do not commute 

Now, to diagonalize at the semi-classical approx- 
imation, it is shown in [23| that it is enough to apply the 
following FW transformation (R being now an operator) 



!7 W (P) = U [r) {P) + X^ 



with 



x (r) 



I 

4fi 



A {r) A (r) 



U (r \P) 



(7) 



(8) 



where we defined the position and momentum pure gauge 



Berry potential A 



0) 

R 



ihU^VpU^ + and A 



(r) 



— i?i{y( r 'VR,£A r ) + , and then to project the transformed 
Hamiltonian JJ^H^U^ + on the positive energy states. 
All expressions in U^(P) are implicitly assumed to be 
symmetrized in P and R and the corrective term X^ 
must be added to restore the unitarity of U^ r \P) which 
is destroyed by the symmetrization. After projection 
on the positive energy subspace needed to perform the 
diagonalization [23j the resulting position and momen- 
tum operators can thus be written r( r ) 



R + A^ } and 



p(r) _ p _|_ w j iere ^- ne ex pii c it computation gives for 
the components 



A 



(r) 
P,k 



% 2 e a ^P a cr {r \v Rk P 7 ) 
-Tic — 



2£('') 2 



+ 0(h 2 ) (9) 



.4 



(r) 



he 2 



2£(r-)2 



(r) 



0(h 2 ) 



(10) 



with £?W the same as E^ above but now R is an oper- 
ator. Performing our diagonalization process leads us to 
the following expression for the positive energy operator 

^(p.rjp.fTW 



c\\v 
ih 



4eM 



V fil (a«.P),V pl (a (r) .P) 



P 



+(V Rl e^)A% + (y plE W)A<$ + 0(H 2 ) (11) 
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Some computations allow us to rewrite the r.h.s. of Eq. 
II II in a more familiar form. Define first the 7-component 
of the vector 

7tp.cr (r) 



e e/ 3 7 rf^(r) and the helicity 



Introducing also 



,(>-) 



< 5 * + ^ — ) ' I l'i 



AWr,(r).p 



P 



an explicit computation shows that the semiclassical en- 
ergy reads 



£ ~M = £ {r) 



pT hB.a^ 

4 p + 2e( r ) 



(A? 



xp).B 



r(r) 



where we have introduced a field B~, = — ^P/iT a ^ s e. 
with T af3S 



h la dih k P - h ll3 dih ka )+h la T 



a/37 



the usual torsion for a static metric (where only space 
indices in the summations give non zero contributions). 

Interestingly, this semi-classical Hamiltonian presents 
formally the same form as the one of a Dirac particle in 
a true external magnetic field [4] 23]. The term B.er is 
responsible for the Stern-Gerlach effect, and the opera- 
tor L = (A^xp) is the intrinsic angular momentum of 
semiclassical particles. The same contribution appears 
also in the context of the semiclassical behavior of Bloch 
electrons (spinless) in an external magnetic field [Bl|(25l| 
where it corresponds to a magnetization term. Because 
of this analogy and since T a @ s is directly related to the 
torsion of space through T a @ s 
a magnetotorsion field. 

However, this form for the energy presents the default 
to involve the spin rather than the helicity, this last quan- 
tity being more fundamental for a photon. Actually one 
can use the property Ap/2p = ha/2 — (A^xp) to rewrite 
the energy as 



h 5 k h ia h^Th we call B 



,(r) 



am p.r aw b. p 



p 



2 £ M 



P 



(12) 



We can now build the Hamiltonian as the sum of the 
two Hamiltonians for one-half massless spinning particle 
e = (e^ + i^)/2. By Taylor expanding the expression 
of we see that at the leading order in % the sum does 
recombine to give after reintroducing the 500 dependence 



£ + A P T o , A .goo B -P 
4 p 2e p 



(13) 



where e 



4 p 



with A = (AW + A (2) )/2 the photon helicity. In the same 
manner, the Berry connections for the photon become 



(i) 



A ( R 2) and A P = A 



(i) 



(2) 



This allow us 



to write the dynamical operators at leading order in h as 



ihdp + he 



P x £ 



ac 2 (^)v R p 



(14) 
(15) 



where in Eq dl4l and [15] e 2 can be approximated by 
c 2 Pi9 Vj 'gooPj at the order considered. 

The semi-classical Hamiltonian EqQ~3] is one of the 
main results of this paper. It contains, in addition to 
the energy term e, new contributions due to the Berry 
connections. Indeed, Eq[l3] shows that the helicity cou- 
ples to the gravitational field through the magnetotorsion 
field B which is non-zero for a space with torsion. As a 
consequence, a hypothetical torsion of space may be re- 
vealed through thepresence of this coupling. Note that, 
in agrement with [24| , this Hamiltonian does not contain 
the spin-gravity coupling term S.Vgoo predicted in [26|. 

From Eqs ll4l and ll5l we deduce the new (non-canonical) 
commutations rules 



[r\r 3 ] = ihQ% 



[p i y] 

[p\r 3 ] = -ihg ij +ih&% r 



(16) 
(17) 
(18) 



where 9£L 
either r or p. 
leading order 



d fi A„ 



-[A0 A„j] where £, 77. mean 
An explicit computation shows that at 



-d v iA Cl 



en = -he 



(£-P)P7_a 



2e 4 



a "73 



■ _ 4 (S.p)p 7 



(19) 



From the additional commutation relations between the 
helicity and the dynamical operators [r^ , A] = [pi , A] =0 
we deduce the semiclassical equations of motion 



d-e, 



V p £- 



p x 8 rr 



(1 - 6 pr ) V r e + r x 9j: 



(20) 



To complete the dynamical description of the photon no- 
tice that at the leading order the helicity A is not changed 
by the unitary transformation which diagonalizes the 
Hamiltonian so that it can be written A = hp.Y^/p. Af- 
ter a short computation one can check that the helicity 
is always conserved 



di V P 







(21) 



for an arbitrary static gravitational field independently 
of the existence of a torsion of space. 

EqsJ^D] are the new semiclassical equations of motion 
for a photon in a static gravitational field. They describe 
the ray trajectory of light in the first approximation of 
geometrical optics (GO). (In GO it is common to work 
with dimcnsionless momentum operator p ^fcp^k with 
fco = uj I c instead of the momentum [ill]). For zero Berry 
curvatures we obtain the well known zero order approxi- 
mation of GO and photons follow the null geodesic. The 
velocity equation contains the by now well known anoma- 
lous contribution p x Q rr which is at the origin of the in- 
trinsic spin Hall effect (or Magnus effect) of the photon in 
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an isotropic inhomogeneous medium of refractive index 
n(r) d, dUdl, ■ Indeed, this term causes an additional 
displacement of photons of distinct helicity in opposite 
directions orthogonally to the ray. Consequently, we pre- 
dict gravitational birefringence since photons with dis- 
tinct helicities follow different geodesies. In comparison 
to the usual velocity f = V p e ~ c, the anomalous veloc- 
ity term is obviously small, its order v]_ ~ cAV r j 
being proportional to the wave length A. 

The momentum equation presents the dual expres- 
sion r x pp of the anomalous velocity which is a kind 
of Lorentz force which being of order Ti does not influ- 
ence the velocity equation at order h. Note that similar 
equations of motion with dual contributions p x O rr and 
f x PP were predicted for the wave-packets dynamics of 
spinless electrons in crystals subject to small perturba- 
tions . The complicated Eq. |2T)1 simplifies greatly for 
a symmetric gravitational field as shown below. 



IV. SYMMETRIC GRAVITATIONAL FIELD 

As a simple application, consider the symmetric case 
9oo9 1 '' = F 2 (R). A typical example of such a metric 
is the Schwarzschild space-time in isotropic coordinates. 
For a symmetric metric one has B.p = F = and the 
semiclassical energy Eq. Q2] reduces to 



(22) 



e = c(pF(r)+F(r)p)/2 

with the dynamical variables r = R + h P p .p , 
the following commutation relations [r 2 ,^] 
—ih\e t: ' k pk/p 3 , [p 1 ^] — 0, [r l ,pi] — %hg % K As a conse- 
quence, we derive the following equations of motion 



p = P, and 

= im% = 



gravitational birefringence predicted here is simply due 
to the Magnus effect as a consequence of the photon spin- 
orbit interaction. In particular this effect does not need a 
coupling between the electromagnetic field and a torsion 
term as proposed in [28[ . 

We now apply the equations of motion to compute the 
deflection of polarized light by a star's gravitational field. 
A polarization independent result is expected by the Ein- 
stein's theory of gravitation which does not consider the 
anomalous velocity. With the Schwarzschild metric one 
has F(R) = 1 - and for the equations of motion 

we get 

P = -2GM^p 

p 2GM r x p 
f = —F + A — ^ 



P 



p- 



We can therefore easily compute the angle of deflection 
A(f> between in-and outgoing polarized rays 



. , AGM I XX 
Ac/)=— (I--; 



c 2 r 



where ro is the smallest distance of the light ray to the 
central source of gravitation, M the mass of the star, A 
the wave length of the photon and A the helicity. 

We observe that the deviation is helicity dependent as 
a consequence of the anomalous velocity, but that the 
effect is very small being of order X/tq and certainly un- 
observable for a star like the Sun. 



CONCLUSION 



= V p e 
= -V r e 



px0 r 



(23) 



In the symmetric case the equations of motion become 
simpler than in the general case, but the gravitational 
birefringence is still present. These equations were al- 
ready postulated (but not derived) in [3] to explain the 
Magnus effect (the different deviation of light of distinct 
polarization in an inhomogeneous medium of refractive 
index n(r)) observed despite its smallness in inhomo- 
geneous isotropic optical fibers [27| and also discussed 
theoretically in less general contexts and with different 
approaches in several other papers [ill [12] 0- This 
case fits within our formalism since a gravitational field 
can be seen as an isotropic medium related to the met- 
ric through the relation = 5* J 'n — 1 (r). Therefore the 



In this paper, we have diagonalized at the first order 
in h the photon Hamiltonian in a static gravitational 
field. This diagonal Hamiltonian displays a new inter- 
action between the helicity and a torsion field. As a con- 
sequence the torsion of space, if any, could in principle 
be determined through this coupling. However, even in 
the absence of torsion, we found two new semiclassical 
equations of motion including Berry phase contributions 
for both dynamical variables, predicting that the photon 
does not follow the null geodesic due to its spinning na- 
ture. The reason is an anomalous velocity, responsible 
for the gravitational birefringence. This last result is in 
agreement with the modern point of view about the spin- 
ning particles evolution. Our results are not restricted 
to the gravitational field but also apply to systems with 
anisotropic refractive indices. 
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